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Abstract. This note is devoted to construct a rough path above a multidimensional 
fractional Brownian motion B with any Hurst parameter H e (0,1), by means of its 
representation as a Volterra Gaussian process. This approach yields some algebraic and 
computational simplifications with respect to [18], where the construction of a rough 
path over B was first introduced. 



1. Introduction 



Rough paths analysis is a theory introduced by Terry Lyons in the pioneering paper [H] 
which aims to solve differential equations driven by functions with finite variation with 
p > 1, or by Holder continuous functions of order 7 e (0, 1). One possible shortcut to the 
. rough path theory is the following summary (see [Tj [HI [H [12] for a complete construction) . 

Given a 7-H61der d-dimensional process X = {X{1), . . . , X{d)) defined on an arbitrary 
interval [0, T] , assume that one can define some iterated integrals of the form 



X^^{ii,. . . = / dX^^{ii)dXu2{i2) ■ ■ -dX^^lQ, 

J S<Ul <-- <u„<t 



^ . for < s < t < T, 77, < [I/7J and Zi, . . . ,z„ G {1, • • • ,d}. As long as X is a nonsmooth 

[ function, the integral above cannot be defined rigorously in the Riemann sense (and not 

even in the Riemann-Stieltjes sense if 7 < 1/2). However, it is reasonable to assume that 
some elements X" can be constructed, sharing the following three properties with usual 
. iterated integrals (here and in the sequel, we denote by Sk,T = {(^i, . . . ,Un) : < ui < 

O I ■ ■ ■ < Un < T} the k^^ order simplex on [0, T]): 

O ■ (1) Regularity: Each component of X" is n7-H61der continuous (in the sense of the 

Holder norm introduced in Q) for all n < [I/7J, and Xj^ = Xt — Xs- 
(2) Multiplicativity: Letting ((5X")s„j := X" — X"„ — X"^ for {s,u,t) G S^^t, one 



X 

^ ■ requires 



n-l 

(5X")s„j(2i, . . . ,in) = (^1) • • • ) "^nJX^i ■ ■ ■ , '"nj 



ni=l 

(3) Geometricity: For any n,m such that n + m < [l/7jand {s,t) G S2,t, we have: 
X^,(zi,...,z„)X-(ji,...,j^)= Yl X^+-(A;i,...,fc„+^), (2) 
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where, for two tuples ^{ij) stands for the set of permutations of the indices 
contained in («, j), and Sh(j,j) is a subset of S(j j) defined by: 

Sh(z, j) = {a G a does not change the orderings of z and j} . 

We shall call the family {X"; n < [I/7J} a rough path over X (it is also referred to as 
the signature of X in [7j). 

Once a rough path over X is defined, the theory described in [71 El [12] can be seen as 
a procedure which allows us to construct, starting from the family {X"; n < [I/7J}, the 
complete stack {X"; n > 1}. Furthermore, with the rough path over X in hand, one can 
also define rigorously and solve differential equations driven by X. 

The above general framework leads thus naturally to the question of a rough path 
construction for standard stochastic processes. The first example one may have in mind 
concerning this issue is arguably the case of a d-dimensional fractional Brownian motion 
(fBm) B = (5(1), . . . , B{d)) with Hurst parameter H G (0, 1). This is a Gaussian process 
with zero mean whose components are independent and with covariance function given 

by 

E {B,{i)B,{i)) = 1 (t^^ + s'"" -\t- sp^) , s,tE M+. 

For H = ^ this is just the usual Brownian motion. For any H G (0, 1), the variance of 
the increments of B is then given by 

E [{Bt{^ - Bsi^f] ={t- (s, t) G S2,T, i = l,...,d, 

and this implies that almost surely the trajectories of the fBm are 7-H61der continuous 
for any 7 < if, which justifies the fact that the fBm is the canonical example for a rough 
path construction. 

The first successful rough path analysis for B has been implemented in \5\ by means of 
a linearization of the fBm path, and it leads to the construction of a family {B-*^, B^, B^} 
satisfying ^ and for any H > 1/4 (see also [6] for a generalized framework). Some 
other constructions can be found in [131 US] by means of stochastic analysis methods, and 
in [T7j thanks to complex analysis tools. In all those cases, the barrier H > 1/4 remains, 
and it has long been believed that this was a natural boundary, in terms of regularity, for 
an accurate rough path construction. 

Let us describe now several recent attempts to go beyond the threshold H = 1/4. 
First, the complex analysis methods used in [16] allowed the authors to build a rough 
path above a process F called analytic fBm, which is a complex-valued process whose 
real and imaginary parts are fBm, for any value of if G (0, 1). It should be mentioned 
however that JJF and Q^F are not independent, and thus the arguments in |l6j cannot be 
extrapolated to the real- valued fBm. Then, a series of brilliant ideas developed in [IBl [H] 
lead to the rough path construction in the real-valued case. We will try now to summarize 
briefly, in very vague terms, this series of ideas (see Section [3] for a more detailed didactic 
explanation): 

(i) Consider a smooth approximation i?^ of the fBm B, and the corresponding approxi- 
mation B"'^ of B". Clearly B"'"^ satisfles relation ([H), but may diverge as e — ^ whenever 
H < 1/4. Then, one can decompose B"'"^ as B"'^ = A"/ + C"'"^, where C"'^ is the 
increment of a function /, namely C^f = ft — fs, and A"'^ is obtained as a boundary 
term in the integrals defining B"'"^. As explained in Section [3l a typical example of such 
a decomposition is given (for n = 2) by A^^^ = —B^ ® 5B^^ and C'^f = B^ ® dB^, 



and in this case ft = /J (g) dB^. Then it can be easily checked, thanks to the relation 
C^f = f^- f^^ that C"f - C^;^ - C"f = for any (s, u, t) e S^^t- This means that C"'^ 
does not affect the multiplicative property ([T]) of B"'^. On the other hand, the boundary 
term A^f is usually easily seen to be convergent as e — > to some limit A". Then, 
the limit A" should fulfill the desired multiplicative property, but it does not exhibit the 
desired Holder regularity k'j for any j < H. It should also be noticed that A"/ is not 
the only natural increment sharing the multiplicative property with B"'"^. We refer to 
Section [3] for further details, but let us mention that another possibility for n = 2 is the 
boundary term SX^^ ® Xf, which is easily seen to satisfy relation (P). 

(a) The essential point in Unterberger's method is then the following. Consider a series 
representation of the fBm B = J2k'^k- We can carry out the above program for each 
component B^ by choosing a particular boundary term. Then, it turns out that there is 
a choice of the boundary term for each component B^ such that their sum satisfies the 
desired Holder and multiplicative properties. This idea has been successfully implemented 
in [TSl [To] using an entire series representation, providing a construction of a rough path 
associated to B. However, this construction is rather long and intricate, first because 
the entire series representation is obtained by a deterministic change of variable involving 
the Cayley's transform and second because the changes in the order of integration in 
the multiple integrals are coded by admissible cuts in some trees associated to multiple 
integrals. This language, well-known by theoretical Physicits (see e.g. [4], or [9] in the 
rough paths context), may sound however difficult to the noninitiated reader. Another 
construction involving the harmonizable representation of fBm is outlined in [20], but this 
method does not avoid the use of tree-based expansions. 

The purpose of the current paper is to take up the program initiated in |18| . and 
construct a rough path over i? in a rather simple way, using the stochastic integral repre- 
sentation of the fBm as a Volterra Gaussian process. We know that (see [HI Proposition 
5.1.3] for a justification) for H < 1/2, each component B{i) of B can be written as 

Bt{i) = [ K{t,u)dWti^, t>0, (3) 
Jr 

where W = (W{1), . . . , W{d)) is a rf-dimensional Wiener process, and where the Volterra- 
type kernel K is defined on M+ x by 



K{t,u) 



ch 



u 



H 



-H 



.H 



2 [V 



\H- 



2 dv 



40<u<i}' 



(4) 



with a strictly positive constant c^, whose exact value is irrelevant for our purposes. 
Then, we show that the simple trick described at point (ii) above can be applied in a 
straightforward way using the Volterra representation, leading to a simple general formula 
for the multiple integrals B". To be more specific, let us describe the main result of this 
paper. 

Theorem 1.1. Let B be a d-dimensional fractional Brownian motion with Hurst parame- 
ter H e (0, 1/2), admitting representation ([^j- For 2 < n < [1/H\, any tuple {ii, . . . ,in) 



of elements of {1, . . . , d}, 1 < j < n and (s, t) G S2^t, set 
B"'^(2i,...,i„) 

„ j-l n 

= (-irW YlK{s,ui)[K{t,u,)-K{s,u,)] 11 K{t,ui) dWu,{ii)---dW^„{i^), (5) 

^^3 1 = 1 l=j + l 

where the kernel K is given by and A" is the subset of defined by 

A'J = {{ui, . . . ,Un) & [0, t]"; Uj = min(Mi, . . . , mi > ■ ■ ■ > Uj_i, and Uj+i < ■■ ■ < m„} . 

Notice that the multiple stochastic integral in ^ is understood in the Stratonovich sense, 
and is well-defined as a random variable as long as n < \1/ H\ . Set also B^j(z) = 

Bt{i)-Bs{i), andfor2<n< ll/H\, 

n 

i=i 

Then the family {B"; 1 < n < [1/H\} defines a rough path over B, in the sense that B" 
is almost surely wy -Holder continuous for any 7 < iJ, and that it satisfies relations (QP 
and 

As announced above, formula ([5]) defines in a compact and simple way the (substitute 
to) iterated integrals of B with respect to itself. Furthermore, this formula also yields a 
reasonably short way to estimate the moments of B" , and thus its Holder regularity. It 
should be mentioned however that our construction is not as general as the one proposed 
in [I9], though it can be extended to a general class of Gaussian Volterra processes (see 
Remark EH). 

In the case 1/4 < H < 1/2, let us also say a word about the relationship between the 
processes B" we have produced and the ones constructed in the aforementioned references 
[SI [131 [E] , which shall be denoted by B"'P (where p stands for pathwise). First of all, let us 
recall that only n = 2 has to be considered for 1/3 < H < 1/2, while n = 2, 3 corresponds 
to the rougher case 1/4 < H < 1/3. Then an easy comparison can be established for 
H = 1/2. Indeed, a slight extension of our construction also allows to define B^ for the 
usual Brownian motion, and it is readily checked in this case that B^ coincides with the 
usual Stratonovich Levy area. Whenever 1/4 < if < 1/2, the situation is less clear: one 
the one hand, we know that (5B" = (5B"'P, and it can be seen from this relation that 
B" and B"'P only differ by the increment of a function On the other hand, it is still 
an open question for us to clearly identify the function We are thus left with two 
candidates B" and B"'P for the rough path construction, B"'P having the advantage of 
being produced as a limit taken on some smooth approximations of the fBm path. 

Here is how our article is divided: some preliminary results, including algebraic in- 
tegration vocabulary, some estimates on the kernel K and Ito-Stratonovich corrections, 
are given in Section [2l Then the basic ideas of the construction are implemented in 
Section [3] on second order iterated integrals. This section is thus intended as a didactic 
introduction to the construction, and could be enough for a first quick glimpse at the 
topic. Then we give all the details concerning the general iterated integral definition and 
prove Theorem 11.11 in Section [H 



2. Preliminaries 



This section is first devoted to recall some notational conventions for a special subset 
(called set of increments) of functions of several variables. These conventions are taken 
from the algebraic integration theory as explained in [HI \10\. We will then recall some 
basic estimates on iterated Stratonovich integrals with respect to the Wiener process, 
which turn out to be useful for the remainder of the article. 

2.1. Some algebraic integration vocabulary. The current section is not intended 
as an introduction to algebraic integration, which would be useless for our purposes. 
However, we shall use in the sequel some notation taken from this method of rough paths 
analysis, and we shall proceed to recall them now. 

The algebraic integration setting is based on the notion of increment, together with 
an elementary operator S acting on them. The notion of increment can be introduced in 
the following way: for an arbitrary real number T > 0, a vector space V, and an integer 
A; > 1, we denote by Sk,T the k*^ order simplex on [0, T], and by CkiV) the set of continuous 
functions g : Sk,T V such that gti - tk = whenever ti = tj+i for some i < k — 1. Such 
a function will be called a {k — 1) -increment, and we will set C^iV) = Uk>iCkiV) . The 
operator 6 alluded to above can be seen as an operator acting on /c-increments, and is 
defined as follows on CkiV): 

k+l 

6 : CkiV) ^ Ck+.iV), i6g)t,...t,^, = J2i-^y9t-r-u-t,^,, (6) 

i=l 

where tj means that this particular argument is omitted. Then a fundamental property 
of 6, which is easily verified, is that 66 = 0, where 66 is considered as an operator from 
CkiV) to Ck+2iV). We will denote ZCkiV) = CkiV) n Kev6. 

Some simple examples of actions of 6, which will be the ones we will really use through- 
out the paper, are obtained by letting g E Ci and h G C2. Then, for any is,u,t) G S^^t, 
we have 

iSg)st = gt- Qs, and i6h)sut = h^t - h^u - Kt, (7) 

and in this particular case, it can be trivially checked that for any g E Ci, one has 66g = 0. 
Conversely, any h G ZC2 can be written as h = 6g for an element g E Ci. In the sequel of 
the paper, we shall write for two elements h^yh"^ G C2 

h} h^, iff = h^ + z with z G ZC2. (8) 

Otherwise stated, iff 6h^ = 6h^. 

Notice that our future discussions will rely on some analytical assumptions made on 
elements of CkiV). Suppose V is equipped with a norm | ■ |. We measure the size of the 
increments by Holder norms defined in the following way: for g G C2iV) let 

\\9h= sup J^'*^ , and C^(y) = {^f G C2(V'); ||5f||^ < cx)} . (9) 
{s,t)&S2,T F - 

With this notation, we also set Cf(V^) = {/ G CiiV); < 00} (notice that the sup 

norm of / is not taken into account in this definition). In the same way, for h G C^iV), 



set 



7,P 

(s,M,t)e53,T 



sup -. -, (10) 



(11) 



where the last infimum is taken over all sequences {hi G C^iV)} such that h = /ij and 
for all choices of the numbers pi G (0, p). Then ||-||^ is easily seen to be a norm on Cs^V), 
and we set 

C^{V) := {heCsiV); ||/it<oo}. 
In order to avoid ambiguities, we shall denote by Af[f; Cj{V)] the //-Holder norm (or 
semi-norm) on the space CjiV), for j = 1, 2, 3. 

In order to define our increments B" for all s, t almost surely, we will also need the 
following definition, which has already been introduced in [161 Section 1.3]: for j > 1 and 
/? > 0, let C™'^(y), where m stands for multiparametric, be the subspace of Cj{V) induced 
by the semi-norm: 

Af[h- Cf^^iV)] = sup |^!^!i±£^^fi^zi!firiif^. e G [0, 1], .1, . . . , G [0,T]| . (12) 

Let us also mention that when = M, we will simply denote the spaces Cj{V),Cj(y) and 
Cj^''^(y) by Cj,Cj,Cf''^ respectively. 

Finally the lemma below, borrowed from [HI Lemma 4], will be an essential tool for the 
analysis of Holder type regularity of our increments: 

Lemma 2.1. Let n > and p > I. Let R G C2(M'), with 6R G C3 (M') in the sense given 

by (dy. // 

■ dudv < 00, 



' S2,T 





Ruv 


|2p 


\u 


— v\ 


2ftp+4 



then R G C2 (ffi ). In particular, there exists a constant Ck,p,i > 0, such that 
Af[R]C^{R^)] 



1 

(f I/? \ 



2.2. Analytic bounds on the fractional Brownian kernel. We gather in this section 
some technical bounds on the kernel K involved in the Volterra representation of B, for 
which we use the following convention (valid until the end of the article): for two positive 
quantities a and b, we write a < 6 whenever there exists a universal constant C such that 
a < Cb. 

First, a classical bound on K is the following: 
Lemma 2.2. Let K be the fBm kernel defined by Then for any < u <t, one has 

K{t,u)<{t-u)''-l+u''-K and d,K{t,u) < {^^y~" {t - uf-l . (13) 

The following simple integral estimate on K also turns out to be useful: 
Lemma 2.3. LetO<v <t<T. Then jl K'^{t, w) dw < {t - vf^. 



Proof. Invoking the bound (fT3l) on K, we have 



K^{t,w)dw< / \{t - w)^-^ + w^-^ 



dw 



< 



2H~1 



dw 



f dw<{t- vf" - [t^^ - v^""] . 



Furthermore, since a" — 6° < (a — fe)" for any < 6 < a and a G (0, 1), we end up with 
K'^{t, w) dw < (t — f )^^, which is our claim. 

□ 

We shall also use a slightly more elaborated result on K: 
Lemma 2.4. Let < s < t < T , assume H < 1/2, and consider the quantity 

1st = j [Kit, ui) - K{s, ui)f i^j K'^it, U2) du^ dui, 

where we recall that we have used the convention K{t,u) = K{t,u)l[Q^t\{u) . Then \Ist\ ^ 



t- s 



\4H 



Proof. According to the fact that K(t, u) = whenever u >t, we obtain the expression 
Ist = [K{t,ui)-K{s,ui)f (^j^ K\t,U2)du2^dui 

+ ^ K\t, Ml) (^j K\t, U2) dw^ dui := II^ + 4- 

Let us bound now the first of those terms: thanks to Lemma [273l one can write J^^ K'^{t, U2) 
du2 ^ (t ~ UiY^ . Moreover, for < -u < s the bound (fT3l) on dtK{t,u) yields 



\K{t,u) - K{s,u) 



dvK{v, u)dv 



< 



(14) 



and thus, putting these two estimates together, we obtain: 

2 







{t - uf^du. 



Performing the changes of variable v = s — u and y = v/(t — s), we end up with 



Ilt<{t-s 



s/(t-s) 



[i + y) 



H- 



.H- 



:i + yy'' dy 



Furthermore, it is easily checked that (1 + y) 



2 — y 2 



1 + y) dy is a conver- 



gent integral whenever H < 1/2, which gives the desired bound for The term I^^ is in 
fact easier to handle, and we leave those details to the reader for the sake of conciseness. 
Then, the estimates on Jjj and J^j yield our claim. 

□ 

Finally, the following related integral bound also turns out to be an important estimate 
for the analysis of n**^ order iterated integrals: 



Lemma 2.5. Suppose that 2kH < 1. For A > 0, set 

-A 



y) 



H- 



H- 



y 



2{k-l)H 



dy. 



Then sup^->q/5a < oo. 

Proof. We can write /3a < «a + 7a, with 



"A 



7a 



y 



dy 



y 



^"^ - (1 + yf-'A (A - 



One can check easily, as in the proof of Lemma [2?H that a a is finite as long as 2kH < 1. 
On the other hand, an obvious change of variables yields 



7A = A 



2kH 



hAiy)il-yr~-^y 



H~\ 2{k~l)H 



dy, 



(15) 



where the (positive) function Ha is defined on IR+ hj hA{y) = y^ ^ ~ + Z/)^ ^- We 
now use two elementary estimates: 



hA{y) < 

and we obtain 

hA{y) = hAiyr"" hAiyy-'''' 

< 

where CH,k = G 



H 



y ^ 
A ' 



and hA{y) < 2/^ 2. 



H 



A 



2kH 



y 



(H-\)(l-2kH) _ '^H,k y 



(l-2k)H-\ 



J^2kH 



H) . Plugging this bound into f fTSl) . we get 



lA < CH,k / (1 - y)''-'^y-''--- 

Jo 



dy. 



This last integral being finite, our claim is now proved. 



□ 



Remark 2.6. The reader can check that the only properties of the kernel K used in the 
sequel are 

K{t, u)<{t- u)^-^ + u^~^ , and dtK{t, u) < {t - u)"'! 

for some HE (0, 1/2). Our rough path construction is thus valid for any Gaussian Volterra 
process defined by a kernel satisfying the above estimates. 

2.3. Contraction of Stratonovich iterated integrals. An important tool in our anal- 
ysis of iterated integrals will be a general formula of Ito-Stratonovich corrections for iter- 
ated integrals. This kind of result has already been obtained in the literature, and for our 
purposes, it will be enough to use a particular case of [2 Proposition 1], recalled here for 
further use. Note that we need an additional notation for this intermediate result: we set 
dY for the Stratonovich type differential with respect to a process Y, while the Ito type 
differential is denoted by dY. 



Proposition 2.7. Let Y = {Y[l), . . . , Y{n)) be a n- dimensional martingale of Gaussian 
type, defined on an interval of the form Yu{i) = ipv{j) dW^i^ij) for a family of 

L^([s, t]) functions {ip{l), ■ ■ ■ , ip{n)), a set of indices (zi, . . . , z„) belonging to {1, . . . , d}^ , 
and where we recall that {W{1), . . . , W{d)) is a d-dimensional Wiener process. Then the 
following decomposition holds true: 



Uj = n 



/I 
dYu^{ii) ■ ■ ■ dYu^in) = XI XI "^^^ 

_<Mi<-<«„<t k=\n/2\ ueDk 

In the above formula, the sets are subsets of {1,2}^ given by 

Dt =<iy= (ni,...,nfc); X' 
and the Ito-type multiple integrals Jst{^) are defined as follows: 

JM= I dz^,{i)---dz^^{k), 

J S<Ul<-<Uf;<t 

where, setting J2'i=i^i = ''^(j); '"^s have 

Z{j) = tf nj = l, 

and 

Zu{j) = (^j^ MmU) - l)^/'„(m(j)) dv^ l(i^(,.)_i=i„o)) if % = 2. 

The previous Ito-Stratonovich decomposition allows us to bound the second order mo- 
ment of iterated Stratonovich integrals in the following way: 

Lemma 2.8. Let ip G Consider the Stratonovich iterated integral 

„ n 

4"M= / ll^{Ui)dWu,{zi)---dWu^{tn). 

J s<m<--<u„<t 



Then, 



^[i:,{^r]<c(^J\{ufdv}j , (16) 

where the constant C depends on n and the multiindex (ii, 

Proof. By Proposition 12.71 we can decompose the Stratonovich integral Istiv) iiito a sum 
of Ito integrals: 

n 

and it suffices to consider each Ito integral Jsti.^)- Then we proceed by recurrence with 
respect to k, with the notation of Proposition I2.7[ Suppose first that rik = 1- Then, 

Jst{,J^) = / Jsu{l^')<^{u)duW{in), 



where z/' = (rii, . . . , rik-i). As a consequence, 

Wst{i^?] = / E [Jsuii^T] V{ufdu < sup E [JsuWf] [ ^{ufdu. 

J S S<U<t Js 



ID* I ^ ^ 



On the other hand, if = 2, then Jst{i^) = jl Jsu{i^")'^{uYdu, with v" = {rii, . . . ,71^-2), 
and again 

E[J.t(z/)'] < sup E [Jsu{v"?] { !\{ufdu 

s<u<t \J s 

By recurrence we obtain (|16p . where C = |^X]fe=[n/2j 2^ 

□ 

3. Iterated integrals of order 2 

In this section, we will define the element announced in Theorem ll.li The study of 
this particular case will (hopefully) allow us to introduce many of the technical ingredients 
needed for the general case in a didactic way. 

3.1. Heuristic considerations. Let us first specify what is meant by an iterated integral 
of order 2: according to the definitions contained in the Introduction, we are searching 
for a process {B^^(ii,i2); {s,t) G S2,t, 1 < H,i2 < d} satisfying: 
(i) The regularity condition G €2"^ {R'^''^) . 
(a) The multiplicative property 

5BL(^l, ^2) = BL(^i) Bl{t2) = - B,{t,)] [^.(z^) - 5„(z2)] , (17) 

which should be satisfied almost surely for all {s,u,t) G S-^^t and 1 < ii,i2 < d. 
(in) The geometric relation, which can be read here as: 

B?,(^i,^2) + B^,(22,^i) = Bl,i^^)Bl,i^2), {s,t) G S2,t, 1 < ^1,^2 < d. (18) 

In order to construct this kind of element, let us start with some heuristic considerations, 
similar to the starting point of [TH|: assume for the moment that X is a smooth d- 
dimensional function defined on [0,T]. Then the natural notion of iterated integral of 
order 2 for X is obviously an element X^, defined in the Riemann sense by 

X?t(?i,^2)= / dX^,{t,)dX,,{t2)= [ [X^{ti)-X,{t,)]dXu{t2)- (19) 

J S<Ul<U2<t Js 

We shall now decompose X^ into terms of the form and as explained in the 
Introduction. In our case, this can be done in two ways: first, equation ffTOl) immediately 
yields 

X^t(^i, ^2) = A^f + C^f , with A^/ = -Xs{ii) SX,t{i2), C^f = / X„(«i) dXu{i2), 

J s 

where we have called those quantities A^'^ and C^'^ because they involve increments of 
the second component X{i2) of X. Notice now that C^'^ is the increment of a function 
/ defined as ft = Xu{ii) dXu{i2)- Hence, according to convention ([8]), one can write 

X^(«i,i2) A^'^. By inverting the order of integration in Mi,M2 thanks to Fubini's 
theorem, we also obtain 

"2 /■ ■ \ _ A 2,1 I r^2,l a2,1 _ c-TA /■ \ V /„• \ /^2,1 



^stih,i2) = Aj +Cj, with A^^ = dX^tiii) Xt{i2), C,i = - / Xu{i2) dX^iii), 
and thus X2(ii,i2) =' A^'i. 



Let us go back now to the case of the d-dimensional fBm B. If we wish the iterated 
integral we are constructing to behave in a similar manner as a Riemann type integral, 
then one should also have the relation: 

B^(zi,22) and B\zut2) A^'\ 

with A^^^ = —Bs{ii) 6Bst{i2) and A^^^ = 6Bst{ii) Bt{i2). This means in particular, ac- 
cording to the fact that S\zc2 = 0) that both A^'^ and A^'^ satisfy the multiplicative 
relation (flTl l. as it can be easily checked by direct computations. However, this naive 
decomposition has an important drawback: the increments A^'^ and A^'^ only belong to 
C2, instead of Cg^, for any 7 < if (this point was also stressed in [T8]). 

Our construction diverges from [H] in the way we cope with the regularity problem 
mentioned above. Indeed, we start from the following observation: invoking the represen- 
tation ^ of B, one can write 

A^f = -Bs{ti)6Bst{t2) = - [ K{s,u,)dWu,{t,) f [K{t,U2)-K{s,U2)]dWu2{t2) 

Jr Jr 

= - f K{s,ui)[K{t,U2)-K{s,U2)]dWuA^i)dWu2{i2), 

where we recall that the stochastic differentials dW are defined in the Stratonovich sense. 
In the same way, we get 

A''/= f [K{t,m) - K{s,ui)]K{t,U2)dWu,{ii)dWu2{t2). 

The idea in order to transform A^'^, A^'^ into increments is then to replace the inte- 
grals over above by integrals on the simplex, as mentioned in the Introduction. Namely, 
we set now 

B^/{ii,i2) = I [K{t,u^)-K{s,u{)]K{t,U2)dW^,{ii)dW^,{i2) (20) 

J Ul<U2 

Bff(zi,Z2) = - / K{s,u^)[K{t,U2)-K{s,U2)]dWuA^^)dW^,{i2), (21) 

J U2<Ul 

and notice that these formulas are a particular case of (l5|) for n = 2. We shall see that 
B^'^(2i,i2) and B^'^(ii,22) are elements of Cg^, but they do not satisfy the multiplicative 
and geometric property anymore. However, it is now easily conceived, by some symmetry 
arguments, that the sum of these last two terms do satisfy the desired algebraic properties 
again. Indeed, we set now 

B?,(^i,^2) =B,Y(«i,«2) + B?f(^i,^2), (22) 

and we claim that B is a Cl^{W^''^) increment which fulfills relations ([ITD and ((IHl). The 
remainder of this section is devoted to prove these claims. 

3.2. Properties of the second order increment. It is obviously essential for the 
following developments to check that B^ is a well defined object in L'^iVL). The next 
proposition asserts the existence of B^^ as a random variable for all s, t in the interval 
[0,T]. 

Proposition 3.1. Let H < 1/2, {s,t) G S2,t andB^^ be the matrix valued random variable 
defined by Then Bl{ti,t2) E L'^{Vt- W^''^) and E[|B?J2] < [t - sf^ . 



Proof. Assume first ii ^ 12- We shall focus on the relation E[(Bfj^(zi, Z2))^] ^ {t — s)^^ , 
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the bound on B^^ being obtained in a similar way. Now Stratonovich and Ito type 
integrals coincide when ii 7^ ^2, and according to expression ( 1201) we have 



[K{t, Ui)l[o,t](ni) - K{s, Ui)liQ^s]{ui)Y K^{t, U2)l[o,t]iu2) duidu2, 

Ul<U2 

which is exactly the quantity 1st studied at Lemma 12^41 The desired bound follows from 
Lemma I2.4[ 

Let us now treat the case ^1 = ^2 = h still concentrating our efi^orts on the inequality 
E[(B^j^(i, i))^] < (t — s)^^. In this context, Proposition 12.71 yields the decomposition 
B^,i\t,t) = Mst + Vst, with 

Mst = I [K{t, ui)-K{s,ui)] K{t, U2) dWu, (0 dW^, {%) 

J U\<U2 

Vst = \ I [K{t,ui) - K{s,ui)]K{t,U2) duidu2, 

^ Jui<U2 

where we stress the fact that Vst is a deterministic correction term. It is thus obviously 
enough to obtain the bounds E[M^j] ^ (t — s)^^ and l^J ^ (t — s)^^ separately, the first 
of these bounds being obtained by evaluating 1st in Lemma 12.41 again. As far as Vst is 
concerned, Jensen's inequality allows us to assert 

Vs] < I [K{t, m) - K{S, U^)f K{t, U2f du, dU2 = 1st, 

which trivially finishes the proof. 



□ 

A second technical step in the study of B^ is to prove that this increment belongs to a 
space of the form C™'^(]R'^''^). 

Proposition 3.2. Let B^ be the iterated integral increment defined by Equation ^2^) . 
Then there exists a version ofB^ such that B^ G C™'^(M'^''^) for any (3 < H. 



Proof. Thanks to the multiparametric version of Kolmogorov's criterion [3], it is enough 
to check that 

EOB?H_.,+,-B?,r] <i/^r+ifcr, (23) 

for (s, t) G i52,T, h, k such that s + h,t + k E [0,T], s + h < t + k, and p large enough. 
Furthermore, invoking the fact that B^^ is an element of the second chaos of W, on which 
all norms are equivalent, it is enough to check relation (l23l ) for p = 2. This last 



computation is however very similar to the one contained in the proof of Proposition 13. 1[ 
and we omit it here for the sake of conciseness. 

□ 

We are now equipped with the continuous version of B^ exhibited in the last proposition, 
on which we will work without further mention, and we are now ready to prove the 
algebraic relations satisfied by our second order increment. 

Proposition 3.3. The increment B^ defined by f2E) satisfies relations [T7\) and ( fTgj) . 



Proof. Recall that we are now dealing with a continuous version of B^. In fact, one can 
easily modify the arguments of Proposition 13.21 in order to get a continuous version of the 
pair (B-'^,B^). This means that it is enough to check relations ( flTll and ([T8ll for some 
fixed 0<s<M<t<r. 

Let us then verify ( JTTl) for s,u,t G [0,T] such that s < u < t. It is readily seen, by 
writing the definitions of B^^^(zi, ^2), B^^^(zi, ^2) and B^^(zi,Z2), that 

5Bl^^Si.^2) = I [K{u,ui) - K{s,ui)] [K{t,U2) - K{u,U2)] dWuA^l)dWu2{^2), 

J Ul<U2 

the right hand side of this equality being well defined as a random variable (a fact 
which can be shown similarly to Proposition 13. II) . Along the same lines, we also get 

6Bt]{^u^2) = [ [K{u,ui) - K{s,ui)] [K{t,U2) - K{u,U2)] dWu,{ii)dWu,{i2), 

J Ul>U2 

and thus 

^BL(^i,^2) = 5B^;t(^i,^2) + ^B^;j(2i,i2) 

[K{u, ni) - K{s, ni)] [K{t, U2) - K{u, U2)] dW^,{zi)dWu2{i2) 



Bl{z,)Bl{z2) 



which is relation ( 1171) . Relation ( iTSll is shown thanks to the same kind of elementary 
considerations, and its proof is left to the reader. 

□ 

Finally, let us close this section by giving the proof of the announced regularity result 
on B^. 

Proposition 3.4. The increment B^ is almost surely an element of €2'^ (Uf''''^) , for any 
-1 <H. 

Proof. Consider a fixed Holder exponent ^ < H. The proof of this result is based on 
Lemma EH which can be read here as J\f[B^; C'^^iR'^''^)] <A + D, with 







2p 




— v\ 


47P+4 



2p 



A=\ I ":L+4 ^^^H ' and D=M[SB^;C 



Let us first deal with the term D above: we have seen that B^ satisfies the multiplicative 
property (fTTll . which can be summarized as 6B^ = 6B ® 6B. Furthermore, B G C1(M.'^) 
for any 7 < i/, and thus, for any I < ii,i2 ^ d and 0<s<u<t<T 

|5BL(^i,^2)| = \5BUii)\ \SBUi2)\ < U^[B; C1{W)] \u - s[<\t - u[< . 

In other words, the quantity ||5B^||^,^ defined by (ITOll is almost surely finite, and according 
to definition (fTTl) . we obtain that D is also almost surely finite. 

We will now show that A is finite almost surely when p is large enough, by proving that 
E[A] < 00. Indeed, invoking Jensen's inequality we obtain: 




where we have used the fact that belongs to the second chaos of W^, on which all 
the LP norms are equivalent. On the other hand, Proposition 13.11 gives [|B^^p] < 
\u — f 1^^^, and plugging this inequality into f l24l) . we obtain that E[/l] is finite as long as 
p>l/{H-^). 

□ 

In conclusion, putting together the last two propositions, we have constructed an ele- 
ment which satisfies the properties (i)-(iii) given at the beginning of the section, for 
any if < 1/2. 

4. General case 

The aim of this section is to prove Theorem 11.11 in its full generality. Recall that we 
define our substitute B" to in}^ order integrals in the following way: for 2 < n < \1/H\, 
any tuple (zi, ■ ■ ■ ,in) of elements of {1, . . . , c?}, 1 < j < n and (s, t) G S2,t, set 



B"'^(zi,...,i„) 




where the kernel K is given by ([4]) and Aj is the subset of [s,t]" defined by 

A] = {{ui, . . . ,Un) e [0, t]"; Uj = min(ui, . . . , u„), Ui > ■ ■ ■ > uj^i, and uj+i < ■ ■ ■ < n„} . 

The increment B" is then given by 

n-l 

B:,{^^,...,^r,) = J2KfiH,...,^n). (26) 

i=i 

It is obviously harder to reproduce the heuristic considerations leading to this expression 
than in Section 13.11 Let us just mention that the same kind of changes in the order 
of integration allows us to produce some increments similar to A^'^,A^'^. Then the 
reordering trick yields some terms of the form B"'-'. After observing the form of several 
of these terms, the general expression (l25l) is then intuited in a natural way. 

Notation: in order to write shorter formulas in the computations below, we use the 
following conventions in the sequel, whenever possible: 

(i) A product of kernels of the form 11^=1 ^i'^j^ ^j) '^ill simply be denoted by 11^=1 -^r, , 
meaning that the variable Uj has to be understood according to the position of the kernel 
K in the product. 

(ii) In the same context, we will also set SK^t for a quantity of the form K{t, Uj) — K{s, Uj). 
(in) Furthermore, when all the Tj are equal to the same instant t, we write YYj=i ^{t^ %) = 

(iv) Finally, we will also shorten the notations for the increments of the Wiener process 
W, and simply write dW for YYj=i^^uj{ij)- 

All these conventions allow us, for instance, to summarize formula (l25l) into 

B^fiz,, . . . , z„) = (-1)^-^ / iff (^--1) dW. (27) 



4.1. Moments of the n**^ order integrals. As in Section [T2l an important step of 
our analysis is a control of the second moment of B". This is given in the following 
proposition. 

Proposition 4.1. For n < [^\, let B"^ be defined by (E^. Then for {s,t) E S2,t, we 
have 

for a strictly positive constant C . 

Proof. Thanks to decomposition (1261 ). it suffices to show that for any fixed family of 
indexes zi, . . . , z„ G {1, . . . , rf} and for any 1 < j < n — 1, we have 

^ . 2" 
^st (^1; • • • ; ^n) 



E 



<C{t-s 



,2nH 



Invoking now expression (l25ll for B"'-^ and decomposing the integral over the region Aj 
appearing in the definition of B"'-'(2i, . . . , in) into integrals over the simplex, it suffices to 
show an inequality of the type 

» n 

E [{Qstf] < C{t - s)2"^, with Qst = / SK,t H JC, dW. (28) 

Notice that in the expression above, we made use of the notation introduced at the 
beginning of the current section, and for i = we assume Tj = s or t. We 

concentrate our efforts now in proving (l28l) . 

Let us further decompose Q into Q = + Q^, where 

Qlt= [ Kf-dW, and Ql= I 5Kst\{K^JW. (29) 

J S<Ul<--<Un<t J 0<Ul<--<Un<t,Ul<S 

Notice that in Ql^ we have assumed the Ti = t for all i, since otherwise this term vanishes. 
Moreover, the term Ql^ can be handled using the properties of the multiple Stratonovich 
integrals established in Lemma 12.81 and applying the estimate obtained in Lemma 12.31 



This yields easily the relation B[{Ql^y] < (t - s)^"^. 
Concerning Q^j, one can write Q^^ = X]j=i ^It where 

„ n 

Bi,= / 5K,t\[Kr^dW. 



i=2 



Notice that in the above equation = t if i = j + 1, . . . , n, since we have again Bl^ = 



otherwise. Each term Bl^ can thus be written as the product of two factors: = ClfD 
where for j > 2 



j 

Ci,= I 5K,t\{K^JW, and Di,= l 

J 0<Ul<--<Uj<S j_2 J S<Uj^l<--<Un<t 



dW, 



and for j = 1, Cj^ = 6KstdW and Dl^ is given by the above formula. 

The random variables Cl^ and Dl^ are independent, and E[(i5^j)^] can be bounded easily 
like E[((5sJ^]. Hence we obtain 

E [{Bif] = E [{Ci,r] E [{Di,f] < CE [(CfJ^j {t - sf^^-^^^ . (30) 



In order to bound the second moment of Cf^, we express this factor as a sum of Ito 
integrals by means of Proposition [221 To do this, we give up for a moment our convention 
on products of increments, and we define, for u G [0, s] and I = 2, . . . , j , the processes 



YJI) 



[K{t,v)-K{s,v)]dW,{ii) and Y^{1) = K{n,v)dW,{zi). 



Then, the processes {¥^{1); < u < s} are Gaussian martingales and 



dY^,{l)dY^,i2)---dY^^il). 



0<«i<-<«,<s 



Thus, a direct application of Proposition 12.71 yields 
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k=lj/2\ i/GDj 



Jos(iy), where JnJiy 



dZ^,{l)---dZ^^{k), 



0<ui<-<ui,<s 



for u = (ji,... , jfc). Thus, setting Ya^iJi = m{h), we have Z{h) = Y{imih)) if jh. = 1, and 
= {Y{m{h) - l),F(m(/i)))^ if jh = 2 and im{h)-i = im{h), where (-, ■) designates 
the bracket of two continuous martingales. We are going to estimate E[Jos(z^)^] using a 
recursive argument. This will be done in several steps: 

Step 1: Suppose jk,jk-i, ■ ■ ■ ,ji = 2. Then j = 2k, and we can assume that im = im-i 
for m = 2, 4, . . . , 2k, otherwise Jos(^) = 0. The term Jos(^) is deterministic and it can be 
expressed as follows: 



'Os{i^) = / [K{t, Ml) - K{s, Ml)] K{t2, Ui) 

k 



X 



JJ K{r2h-i, Uh)K{T2h, Uh)dui ■ ■ ■ duk- 



h=2 



As a consequence, owing to (fTS]) and ([HI), we have 



l^osMI<C / ^^^}\[^^^ldu,---duk 



(31) 



where 



\S - Ui 



-1 



Is — Ml) 2 -|- u 



and 
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Moreover, the integral of 11^=2 V^^^^ easily bounded: indeed, we have 



n '^S duk 

k 



< 



Ul<U2< - <Uk<S 



Uh-1 



<c 



ui<U2< --<Uk<s 



du2 ■ ■ ■ dUk 

du2 ■ ■ ■ duk 



<C(s-Mi)2('="i)^, 

with the convention Uk+i = s. Therefore, plugging this inequality into ( !3Ti ) and making 
the change of variables s — ui = v and y = we get 



c 



v^~^ -{t- s + v)^~^ 



^ + (S - v)^ 5 



s/{t-s) 



t-s 



We are now in a position to use Lemma 12751 with A = s/{t — s), and we obtain 

\Mi^)\<c{t-sf''', 

which implies that Jos(^)^ < C(t — s)^''^ , owing to the fact that 2k = j. 
Step 2: Suppose that jk = 1. Then Proposition 12.71 gives 



(32) 



dZ^,{l) ■ ■ ■ dZ^^_^{k - l)i^(r,-, Uk) dW^{'. 



0)1 



0<ui<---<Uk<s 



and 



E [Miy?] 



E(Jo„(i/')') K{Tj,ufdu 



Jo 



u 



^) du, 



with z/' = (ji, . . . ,jk-i)- This relation allows to set an induction procedure, as we shall 
see later. 

Step 3: Suppose that jk, jk-i, ■ ■ ■ , jb+i = 2 and jb = I, where b > 2. We assume that 
^m(fc) = im{h)-i for h = b + 1, . . . , k. Here again, Proposition 12.71 implies 

k 



,Uh)K{Tra(h),Uh)dUi ■ ■ ■ dUk, 

0<ui<-<Ufe<s h=b+l 



and Fubini's theorem yields 



Jos{jy) = / Jout{jy')K{T„i(^h),Ub)G{Ub)dWu^{im(h)), 



with u' = (ji, . . . , Jb-i), and where 



G{ub) = Yi K{T^(^h)-i,Uh)K{T'^''''\uh)dub+i ■ ■ ■ duk- 

As for the previous bound (!32l) we obtain 

Therefore 

Jo 

< C Te [JouA^'r] [is - Ubf"-' + «f (. - UbY^'~'^^ dUb. 
Jo 

Notice that the above inequality includes the inequality obtained in Step 2, which corre- 
sponds to the case b = k. 

Step 4: Suppose that jk, jk-i, ■ ■ ■ , jb+i = 2, jb = I, jb-i, jb-2, ■ ■ ■ , jc+i = 2 and jc = I, 
where 2 < c < b. We assume also that im{h) = im(h)-i for /i = c+1, . . . , , 6— 1, 6+1, . . . , /c. 
By the same arguments as in Step 2 we obtain 

E [Mvf] < C f E [JouSi^r] [(% - u.r"~' + (ub - 

J 0<Wc<«6<S 

X [{s - Ubf''^' + {s - UbY^'-'^^'duJub, 

with u' = (ji, . . . , ic-i)- Replacing u\^~^ by [ub — Uc)^^~^ and integrating with respect 
to Ub yields 

^0 

Step 5: Iteration scheme. Iterating the argument in Step 4, we reduce the size of z/' until 
we obtain a multiindex of length r such that z/' = (1, 2, . . . , 2) or = (2, 2, . . . , 2), with 
jr+i = 1, and we obtain an estimate of the form 

E [Mi^f] <C Te [Jouii^r] [is - uf^-^ + n^^-i] (. - ur^t.,.ndu^ (33) 
Jo 

Suppose first that u' = (1, 2, . . . , 2). Then, 

W) = I [K{t,ui)-K{s,ui)] 

J 0<U\<--<Ur<U 

r 

X '^K{Tm{h)-l,Uh)K{Tm{h),Uh)dWu^{ii)dU2 " --dUr. 
h=2 

and by Fubini's theorem 

Josiy') = r[Kit,u,) - K{s,u,)]Fiu,)dW^,iH), 
Jo 



where 



-^(^1) = / Y\. ^i'^rn{h)-l, Uh)K{Tm{h), Uh)dU2 ■ ■ ■ dUr- 

J Ul<U2<---<Ur<U 



As in the proof of (l32|) we get 

\F{ui)\ <C(u-ui)2(''-i)^. 

Therefore, 

E [Josii^T] <C [(t - -is- u,f~'^]\u - u,Y^^-'~^" du,. (34) 

^0 

Substituting (l34l ) into ( !33l l yields, after integrating in the variable m, 

E [Myf] < C f[{t - - (s - - 

Performing the changes of variables v = s — u and y = v/(t — s), we end up with 

E [Mi^Y] < C{t - sf^"" / [(1 + t)^-^ - i/^^5]2^20-i)/f^^ < _ ^)2j/f^ (35^ 

Jo 

where the last step is obtained thanks to a slight variation of Lemma I2.5[ 
If u' = (2, 2, . . . , 2), then we proceed as in Step 1 and we obtain 

\JouW)\<C f\{s-u,)'''^^-{t-u,)^-^[{u-u,f-'^+uf^^^^ (36) 

Substituting ( l36l) into ( 1331 ). integrating first in the variable u and using the same arguments 
as in Step 1 we obtain also the estimate 

^[Myf] <C{t-sf^''. (37) 

Step 6: Conclusion. Our bounds ( l35l ) and (l37l) on Jos(i^) yields the same kind of estimate 
for the term Cl^. Thus relation (1301 ) gives Bl^ < (t — s)^"^. This estimate can now be 
plugged into the definition (1291) of Q^, then in the definition of Q, which leads to our claim 
( l28ll . The proof is now finished. 

□ 

4.2. Proof of Theorem II. IL Before we prove our main theorem, we need a last elemen- 
tary technical ingredient, which relies on the notational convention given at the beginning 
of the current section. 

Lemma 4.2. For n > 3, j = 2, . . . ,n — 1 and < s < t < T , set 
Recall that for an element M G C2, SM is defined by ([3j. Then 

m=l 

n-j 
m=l 

The relation still holds true for j G {l,""-} and n = 2, with the convention K®^ = 1 and 
= 0. 



Proof. This proof is completely elementary, and included here for the sake of completeness, 
since it uses heavily the notation of Section I2.1[ 

First, if a, b, c are 3 increments in Ci, and if we define G C2 by N^t = Sbgt Ct, then 
a simple application of Definition ([7]) gives 

SNsut = -Sttsu Shut Ct + as 6bsu 5cut- 

Our claim is thus proved by applying this relation to a = K'^^^~^\ b = K , c = K'^^"~^\ 
and observing that [6K'^%t = Ep=i Kf^"''^ SK^t Kf^'^K 

□ 



Proof of Theorem \l.l[ The structure of the proof is the same as in the second order case 
of Section I3.2[ we first prove that (a modification of) B" is almost surely an element of 
C™'^((M'^)®") for any (3 < H, applying Kolmogorov's criterion and with the same kind of 
computations as for Proposition 14. 1[ This allows us to reduce the algebraic relations ^ 
and (El) to the case of some fixed s,u,t. We will first focus on ([1]). 

Step 1: Proof of the mutliplicative property ([2]j. Fix {s,u,t) G S^^t- Recall that B"'-' is 
defined by (|2711 . Therefore, invoking Lemma [121 5B"'-' is given by 

5B^^(zi, . . . , z.) = (-1)^- / V Kf"'-'^ 5Ks^ Kfj-i-m) ^k^^ k^^.-,) 

^^7 m=l 



On the other hand, set Z^ut = Er=i B^„^B^i""\ One can easily check that 

n— 1 ni n~ni 

Zsut = Y.Y.Y.^^u'Kr"' (39) 



ni=l k=l h=l 
n—1 ni n—rii 



1 7—1 JA^,h{ni) 



ni=l k=l h=l 

where Ak,h{ni) is the set defined by 



Ml, . . . ,M„); Uk < Uk+l < ■ ■ ■ < Uk < < ■ ■ ■ < Ml, 



We want to show that (l39l) and ([38]) coincide. 

In order to follow the computations below, it might be useful to keep in mind an 
illustration of the coordinates ordering on a set of the form A^^hi^), for which an example 
is provided at Figure [H (note that the ordering between Um and Um+i is not specified). 

Notice that on the set Ak,h{ni) fl {uk < ?ini+/i} the minimum of the coordinates is Uk, 
and on the set Ak^hi^i) H < u^} the minimum is Um+h- Define 

^/ifc(^i) = ^fc,/i(?T.i) n {wfc < Mni+/i}, and Ali^{ni) = Ak,h{ni) n {un^+h < Uk}. 



+ 

n 



Figure 1. Coordinates ordering on Ak^h{m)- 
Consider now the decomposition Z = + Z"^ , where 

n— 1 ni n—ni „ 
Zl^ =J2J2Y1 (-1)'^' / ^^'^'"'^ "^^^ j^^in,~k+h-l) ^^^^ j^^in-n,-h) 

ni=l k=l h=l -^4,h("i) 

We fix j and we try to compute the contribution of on the set for i = 1,2. 
This contribution will be the sum of the integrals on the set A^n /l^^(ni), for each 
/c = 1, . . . , TT-i, h = 1, . . . ,n — rii and for each rii = 1, . . . ,n — 1. 

Notice first that the intersection A^n /l^,^(ni) is non empty only if k = j, h = 1, and 
Um < Um+i which also implies j < rii. Moreover, in this case we have A^ -^^{ni) fl {um < 
Um+i} = 111 this way we obtain that the contribution of on Aj is 

n-l 

(-i>'-' E 

ni=j 3 

= (-1)^-1 V / Kf^-^HK^^Kf'^-^UK^tKf^-'Uw, (40) 
m=l ^^7 

where we have used the simple change of variables m = rii — j + 1. In the same manner, 
on the set A^n we have k = rii, h = j , which also implies rii < j. Therefore, the 

contribution of Z]^^ on Aj is 

/ Kf'^'^^UKsuKf'-^-"'^ K^^""''^ dW. (41) 

One can now easily verify that the sum of fllOl) and fllTll is equal to the term (l38l) . 

It remains to prove that the contribution of Zgut to the set {UjA^Y is zero. For this, 
observe that {UjA'jY can be split into slices Dk,p,h of the following form: for 1 < < p < 
n — 1, we assume that Uk < Uk-i < ■ ■ ■ < ui and Uk < Uk+i < ■ ■ ■ < Up but Up > Up+i. 
Suppose also that 1 < h < n — p and that Wp+h is the minimum of the coordinates 
Mp+i, . . . Then, for -Dfc,p,h to be a subset of U^^^^ Uk^h Ak,h{ni), we need the further 
condition Up^h < Up^h+i < ■ ■ ■ < Un and Up+h < Up^h-i < ■ ■ ■ < "Up+i. With all these 
constraints in mind, it is easily seen that D^ p h corresponds to two possible choices of set 
^k,h{ni). Indeed, we have 

Dk,p,h = Ak,h{p) = Ak,h+i{p - !)• 

Going back now to the expression ( 1391 ) of Zgut, it is readily checked that the two contri- 
butions, respectively on Ak^h{p) and Ak^h+i{p — 1), yield two terms with opposite sign, 
which cancel out in the sum. 



Step 2: Proof of the geometric property ([^. Fix n, m such that n + m < [I/7J and let 
(s, t) G S2,T- Consider the product 

B^,(2i,...,^„)B^(ji,...,jJ 

where we have used notation (1271) and where we recall that the sets Aj and are defined 
by 

= {u G [0, t]" : Uj < Uj^i < ■ ■ ■ < Ui,Uj < Uj^i < ■ ■ ■ < 

= {v G [0, t]"" : Vh < Vh-i < ■■■ <vi,vh < Vh+i <■■■ < Vm}- 

The product of the two Stratonovich integrals can be expressed as a Stratonovich integral 
on the region A" x A]^ with respect to the diff'erential 

We will make use of the notation z = {u, v), where ^q, = Mq, for a = 1, . . . , n and Za = Va-n 
for a = n + 1, . . . , n + m. Like in Step 1, the region Aj x A^ can be first decomposed 
into the union of the disjoint regions Dj ^ and Ej h , corresponding respectively to the 
additional constraints {uj < Vh} and {uj > Vh} (notice that this decomposition is valid 
up to the set {uj = Vh}, whose contribution to the stochastic integral is null). 

Consider first the case {uj < Vh}- On Dj^h the minimum of all the coordinates 
Then Dj^h can be further decomposed into the disjoint union of the sets 

Dlh,i = {ze [0, : zj < z^^^,^_^ <---<z^,,Zj<Zp,<---< J 

where 

7r(l, . . . , n + m) = (ai, . . . , aj+h-2, j, A, • • • , ^n-j+i+m-h) 

runs over all permutations of the coordinates 1, . . . ,n + m such that 7r(j + h — l) = j and: 

(i) ai, . . . , aj+h-2 is a permutation of the coordinates 1, ... ,j — 1 and n + 1, . . . ,n + h — 1 
that keeps the orderings zi > ■ ■ ■ > zj^i and Zn+i > ■ ■ ■ > Zn+h~i- 

(a) Pi, . . . , Pn-j+i+m-h is a permutation of the coordinates j + . . . ,n and n+h, . . . , n+m 
that keeps the orderings zj+i < ■ ■ ■ < Zn and Zn+h < ■ ■ ■ < Zn+m- 

Moreover, Dj h can be also be decomposed into the disjoint union of the sets 

Dlh,2 = [0, : Zj < Z^.^^_^ <---<Za,,Zj<Zf3,<---< ^/3„_,+,„. J 

^{Zn+h < Zn+h+l}, 

where 

7f(l, . . .,n + m) = («!,..., aj+h-lj, A, • • • , (3n-j+m-h) 

runs over all permutations of the coordinates 1, . . . ,n + m such that 7r(j + h) = j and: 

(i) ai, . . . , aj^h-i is a permutation of the coordinates 1, ... ,j — 1 and n + 1, . . . ,n + h 
that keeps the orderings zi > ■ ■ ■ > Zj_i and Zn+i > ■ ■ ■ > Zn+h- 

(a) Pi, ... , Pn-j+m-h is a permutation of the coordinates j+1, . . . ,n and n+h+1, . . . , n+m 
that keeps the orderings zj+i < ■ ■ ■ < Zn and Zn+h+i < ■ ■ ■ < Zn+m- 



Then, on the set Dj ^ we write 

= irf (^-i)5i^,,i^f ("-^'^irf (''-i^i^f - Kf^-'^6KstKf''-'^Kf''Kf"'-''\ 

and the integral 

can be expressed as the sum /j ^ = 1'^^^ + /J"^, with 

j+h-2 



n— j'+l+m— h 
1=1 



and 



. j+h-i 



^j,h,2 1=1 

n—j+m—h 

n+m ) • 

1=1 

Let us handle first the term I^^^: consider a permutation a of 1, . . . , n + m which maps 
aj_|_/i_2 into 1, . . . , j + /i — 2 and . . . , (3n-j+i+m-h into j + h, . . . ,n + m, with the 
additional condition cr(j) = j + /i — 1. If we make this permutation in the coordinates of 
^Ih obtain 



. j+h-2 



n+m 

X K{t,Zi)dW,,{k^)---dW,^^^{kn+m). 
l=j+h 



where /ci, . . . , /cn+m is a permutation of the indexes «i, . . . , ^n? Ji? ■ ■ ■ preserving the 
orderings of the indexes ii, . . . , i„ and ji, . . . ,jm, and which puts ij in the j + — 1 place, 
and where u, r] are defined by 

u = mm{i >j + h:kie {ji, . . . , j™}}, r] = max{i < j + h - 2 : ki E {ji, . . .,jm}}- 



In the same way, we can consider a permutation a in the coordinates Zi which maps 
aj+h-i into 1, . . . , j + /i - 1 and . . . , Pn-j+m-h into j + /i + 1, . . . n + m, and 



'^(j) = J + ^- If W6 make this permutation in the coordinates of Ij we obtain 

j+h~i 

(-1^+^-1 H K{s,zi)5K,t{zj+h-i) 

n+m 

X II K{t,zi)dW,,ih)---dW,^^JK+J, 
i=j+h+i 

where /ci, . . . , kn+m is a permutation of the indexes ii, . . . , in,ji, ■ ■ ■ ,jm which does not 
change the orderings of the indexes ii, . . . ,in and and ji, . . . and where z/, 77 are now 
defined by 

u = mm{i > j + h+l : ki e {ji, . . . , jm}}, V = max{i < j + h - 1 : ki e {ji, . . . , jm}}- 

When we sum these integrals over all permutations a of the above type, and also over 
j and h, we obtain Xlfcgsh(fj) B"^™'^/;;!, • • • , kn+m), where 

n+m „ P~l 

p=l,fepG{n, .■■,««} '=1 

n+m 

«=p+i 

In a similar manner we could show that the sum of the integrals over Ej^h give rise to 
X)fcesh(ij) B^t^'^'^l^i, kn+m), for /i = 1, . . . , m, where 

n+m „ p— 1 

B^+'"'^(fci,...,A;.+^) = 5^ / {-l)^-'l[K{s,zi)6K,t{z,) 

P=i,fcpe{ji,...j™.}^p z=i 

n+m 

i=p+i 

Taking into account the two contributions B"/'"'^ and B"/'™'^, the proof of the geometric 
property is now easily finished. 

Step 3: Proof of the regularity property. Like in Proposition 13.41 the fact that B" belongs 
to C2 ^ for any 7 < is an easy consequence of the moment estimate of Proposition 14.11 
plus a simple induction procedure. 

Indeed, assume that B**^ G C2^((]R'^)®^) for any k < n — 1. Then Lemma ETT] gives here 
that Ar[B"; C^\m'^''^)] <A + D, with 

_L_ 

Furthermore, since we have seen that B" satisfies the multiplicative property ([T]), then 
D is easily shown to be almost surely finite thanks to our induction hypothesis. Finally, 
the quantity £[^4] can be bounded along the same lines as in Proposition 13.41 except that 
Proposition 14.11 is used instead of Proposition 13. 1[ 

□ 
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